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ABSTRACT 


This project quantifies several aspects of a new 
multiplexed imaging technique proposed by D.S. Davis. The 
novel approach of this technique involves the use of encoding 
masks derived from a basis set of two-dimensional Walsh 
functions. There were two distinct problems addressed by this 
thesis research. First, a study of computer simulated 
diffraction patterns of the photon flux through these encoding 
masks yielded design constraints to be incorporated into a 
prototype system. These constraints were expressed in a 
simple mathematical relation in terms of wavelength, 
diffraction angle, and spatial frequency. A second problem 
addressed the minimum spatial resolution required for pattern 
recognition. The conclusion reached was that the minimum 
number of resolution elements necessary for pattern 
recognition is 64 in each direction. This determination also 
fixed the minimum size of the basis set of two-dimensional 
Walsh functions required for multiplexing, as well as the 


number of pixels required to display the image. 
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I. INTRODUCTION 


A. MOTIVATION 

This thesis is part of a much larger research and 
development effort being conducted by the author's thesis 
advisor {Ref. 1]. The ultimate goal of this larger project is 
to develop and evaluate a new type of instrument for infrared 
imaging and imaging spectroscopy. The key feature of this new 
technique is that it will multiplex images, by means of focal 
plane encoding masks, onto a single detector, or at most a 
pair of detectors. The viability of this method was 
demonstrated in the proof-of-concept thesis research conducted 
by Capt R.H. McKenzie, III, USMC [Ref. 2]. The motives for 
this thesis project have been twofold. First, since the 
masking technique will be employed at long wavelengths, a 
quantitative prediction of the diffraction behavior is needed 
before mask design is undertaken. This project successfully 
addresses that problem. The second motive for this research 
project is a direct consequence of the multiplexing technique 
itself. The spatial resolution of a multiplexed image will be 
inversely proportional to the number of encoding masks 
incorporated into a working instrument [{Ref. 1]. This 
project, therefore, successfully provides initial estimates 


for the number of encoding masks needed to produce 


recognizable multiplexed images from a variety of commonplace 


objects or targets. 


B. BACKGROUND 

Terrestrial objects emit most of their thermal radiation 
at infrared (IR) wavelengths. Many interesting objects are in 
approximate thermal equilibrium with their surroundings, which 
have temperatures in the 275 - 310 K range. Assuming that 
these sources radiate approximately as blackbodies [Ref. 3], 
the bulk of their radiance lies within the 5 - 25 um 
wavelength range. Of course, warmer, man-made objects will 
radiate at correspondingly shorter wavelengths. It is 
desirable to utilize this abundance of radiated energy to 
detect and to identify objects. Infrared images provide 
information about the size, shape, and temperature of the 
radiating sources. Infrared spectra reveal the chemical 
composition and thermodynamic information. 

D.S. Davis, advisor on this thesis project, has been 
involved with state-of-the-art spectroscopic IR 
instrumentation for many years. Historically, his research 
interests have involved the development and use of unique, 
high resolution, fully multiplexed, infrared Fourier transform 
spectrometers. His interest in IR multiplexing techniques 
inspired the concept of multiplexed imaging, which is one of 
the goals of the larger research project of which this thesis 


is a part. [Ref. 1] 


The purpose of infrared imaging is to measure the image's 
irradiance as a function of two independent spatial variables, 
such as conventional Cartesian (x,y) coordinates. 
Spectroscopic imaging introduces a third independent variable, 
such as wavelength or frequency. Traditional IR imaging 
devices and imaging spectrometers typically scan over one or 
more of the independent variables which describe image 
irradiance. There is, however, an inherent inefficiency 
associated with such scanning techniques. This inefficiency 
stems from the fact that, by observing only one small element 
of the image field (and/or one spectral element) at a time, 
the instrument ignores most of the energy from the scene. 
This can have severe consequences in the infrared, because 
available photons from a target are precious, and they 
typically comprise only a small fraction of the entire, 
background-dominated, IR irradiance. The goal of the proposed 
technique is to develop a more efficient generation of 
instruments which will apply the proven advantages of infrared 
Fourier transform spectroscopy (FTS) to infrared imaging and 


infrared spectroscopy. [Ref. 1] 


C. STATEMENT OF THESIS PROBLEM 

This thesis research addresses two fundamental questions 
which must be answered before the design of a prototype 
multiplexed imaging instrument is complete. First, since the 


instrument will operate at comparatively long wavelengths, 


there is concern that diffraction effects produced by the 
encoding masks will cause substantial beam spreading and 
degradation of image resolution. Therefore, the initial 
portion of this thesis project quantifies the diffraction 
behavior of the masks, so that this behavior can be included 
in subsequent instrumental design considerations. Second, the 
image spatial resolution required to produce identifiable 
images of ordinary terrestrial objects must be determined. 
When the proposed encoding scheme operates at_ short 
wavelengths (far removed from the diffraction limits 
considered in the first part of the research), it is 
conceptually capable of extremely high spatial resolution. In 
fact, the resolution will be inversely proportional to the 
number of encoding masks used. Determination of the minimum 
useful number of masks is important so that needless effort 
and resources are not wasted on an over-designed prototype 


instrument. 


II. OVERVIEW OF MULTIPLEXED IMAGING 


This section provides a synopsis of the multiplexed 
imaging technique. Although the ultimate goal of the entire 
project is to develop a multiplexed imaging spectrometer, the 
thesis research at hand considers only multiplexed imaging. 
Therefore, the spectroscopic aspects of the larger project are 


ignored. 


A. MULTIPLEXING CONCEPTS 

The key feature of optical multiplexing, and the one that 
produces its inherent efficiency, is that all of the photons 
from the target are detected all of the time. This is 
achieved in a multiplexed scanning instrument by detecting 
linear combinations of the signals of interest, rather than by 
scanning through sets of the individual signals themselves. 
For instance, consider the operation of a raster scanning 
imager, such as a FLIR. Let @; be the radiative flux through 
the ith pixel of the image, and P, be the instantaneous power 
sensed by the detector. Assume that there are  0no 
inefficiencies or losses associated with the detection 


process. Then, for N pixels which comprise the entire image, 


iP: = >; ’ {i = 1) «cunt ® 
(1) 


While the device observes the ith pixel, it ignores the other 
N - 1 pixels. 
A multiplexed sensor detects a linear combination of the 


pixel fluxes 


N 
Ee ma S> Hi j9; ’ {1 = 1,..-,N} , 
=a 


(2) 


where the u;,; are weighting coefficients. Therefore, one scans 
through N u's rather than N ¢'s. The multiplexing weights u;, 
are not chosen aebiecac ee there are several constraints 
which they must meet (Ref. 1]. First, equation (2) must be 
invertible, so that the individual pixel fluxes ¢; can be 
recovered from the measurements of the P.. Second, the Mi 
need to permit as much radiation as possible to pass to the 
detector, in order to maximize the overall optical efficiency. 
The inherent inefficiency of traditional raster scan devices 


is due to the fact that their 


Hij = $i; ’ 
(3) 
where 6;; 1s the usual Kronecker delta symbol. Third, the u;, 


must be capable of being implemented in a practical optical 
design which is suitable for a working instrument. 

For example, consider the process by which a Fourier 
transform spectrometer (FTS) encodes spectral intensities 
(Ref. 1]. In this instance, the ¢; represent spectral fluxes 
at N different frequencies or wavelengths. The pw,; are 


discrete samples of the familiar orthogonal sinusoidal 


functions of Fourier analysis. The P, are the discrete 
samples of the recorded interferogram. Inversion of equation 
(2) to yield the spectral fluxes is accomplished by applying 


a discrete Fourier transform to the P.. 


B. SEQUENCY THEORY AND MULTIPLEXED IMAGING 

The family of functions chosen to implement the encoding 
scheme for this multiplexed imaging project are the Walsh 
functions from sequency theory [Refs. 1,4] rather than the 
sinusoidal functions associated with FTS and Fourier analysis. 
The Walsh functions are similar to the sinusoids in that they 
constitute complete, orthogonal basis_ sets. In other 
respects, however, they are radically different. The most 
striking differences are that they are discontinuous, and that 
they assume only one of two possible values, + 1. Another 
difference is that the Walsh functions are not necessarily 
periodic. Furthermore, a complete basis set of Walsh 
functions must contain N = 2" members, where m is a 
nonnegative integer. Detailed discussions of these and other 
properties may be found in Beauchamp [Ref. 4]. 

The concept of a Walsh function's sequency is fundamental 
to understanding this thesis research. Its definition is 
straightforward. The domain over which a basis set of Walsh 
functions is defined is partitioned into N subdomains of equal 
length, where N is the number of functions in the complete 


basis set. The Walsh functions then assume constant values, 


either +1 or -1, across each subdomain, according to the rules 
which dictate their generation [Ref. 4}. At the boundaries 
between adjacent subdomains, a function may switch values, 
thereby crossing zero. The total number of such zero 
crossings that occur across the entire domain of the 
function's definition is that function's sequency. For a 
complete Walsh basis set, the corresponding member function's 
sequencies range from zero toN- 1. N is referred to as the 


order of the basis set. 


C. THE NEW IMAGE ENCODING TECHNIQUE 

The new multiplexed image encoding technique will utilize 
Walsh functions to generate optical encoding masks. These 
masks will, in turn, be used to modulate the flux distribution 
in an image, in analogy with equation (2). An illustration of 
how this might be done is provided by the following example. 
Suppose that one desires to determine the radiant flues ?; 
through different contiguous regions of an image field. 
Consider the complete Walsh basis of order 4. Construct a 
one-by-four row matrix for each basis function in the set, 
such that each function's subdomain structure corresponds to 


one entry, or matrix element: 


Sequency Row Matrix 


0 (+1 +1 #+1 ~°«&421) 
1 (-l1 -1 +1 °~»& 42) 
2 (-l +1 #+1 2 ~«-121) 
3 (+41 -1 +1 = &«©-1) 


The reader may easily verify that the rows are mutually 
orthogonal; in fact, so are the columns. Hence, the rows 
could be combined into a single four-by-four orthogonal Walsh 
encoding matrix, such that 


+1 +1 41 +1 


-1 -1 +1 +1 
“Sea eS. eomeetl|| 


+1 -1 +1 -1 
(4) 


The practical optical implementation of this matrix would 
involve construction of four masks, each of which is portioned 
into four regions where boundaries match those of the image 
regions to be multiplexed. A possible example of such a 
portioning is shown in Figure 1. Each mask in Figure 1 
corresponds to one row of matrix (4). Dark regions imply a 
matrix element of +1 and signify that the mask is optically 
transmitting across those pixels. Light regions correspond to 
matrix elements of -1 and denote that those pixels are 
optically reflective. Light that is transmitted is collected 


and sent to one detector; that which is reflected is sent to 


a second detector. The second detector's output is subtracted 
from that of the first detector. As the encoding instrument 
cycles through the four masks, the differenced detector signal 
will be proportional to the P, in equation (2). Collecting 
the four P,'s into a single column vector and multiplying it 


by the inverse of u;; would yield the desired ¢;'s. 
Figure 1 Two-dimensional Sequency -Ordered Walsh Masks for 
Four Pixels 





The actual encoding scheme in a prototype multiplexed 
imaging system will not make explicit use of the two- 
dimensional mask configurations like those in Figure 1. 
Rather, the prototype will utilize two essentially one- 
dimensional mask families, and it will exploit a very powerful 
mathematical feature of the Walsh functions, the Kronecker 
product [Ref. 4]. Via this approach, a two-dimensional mask 
basis may be synthesized from two lower-order one dimensional 
bases. This is again illustrated by example. Consider the 


Walsh basis of order two. Again, write down the corresponding 


matrices, but this time as both row (one-by-two) and column 


(two-by-one) matrices, as shown below. 


(+1 +1) (-1 +i» 





Multiply together the two numbers located on the outer 
row-column intersections, and collect the products in a two- 
by-two groupings, as shown. Note that these groupings, called 
Kronecker products, are exact numerical analogs of the mask 
pixel structures in Figure 1. The new multiplexed imaging 
technique will implement an optical scheme to incorporate this 
property, thereby greatly reducing the number and complexity 
of the encoding mask families required [Ref. 1]. Figure 2 
illustrates a complete 64 pixel "checkerboard" mask basis, 


synthesized from two one-dimensional bases of order eight. 


D. RESTATEMENT OF THESIS GOALS IN WALSH SEQUENCY-THEORETICAL 
TERMS 
As a consequence of the foregoing discussion, the two 
goals of this thesis research expressed in section I.C. may be 
rephrased in sequency-theoretical terms: (1) since the 
prototype system will generate elaborate two-dimensional image 


encoding masks via the Kronecker product, the system designers 
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need to know how these masks will diffract long wavelength 
radiation, as well as encode it. Diffraction effects could 
cause severe beam spreading problems, resulting in substantial 
loss of optical efficiency, particularly at these longer 
wavelengths and/or when using encoding masks of high spatial 
sequency. Therefore, the first goal of this research has been 
to predict these effects quantitatively. (2) The second goal 
is stated more succinctly: the system designers need to know 
how many spatial sequencies are sufficient to encode images 
which remain recognizable once decoded. It would be both 
expensive and more optically complex to build a prototype 
instrument with a far larger Walsh basis set (i.e., containing 


bases of higher order) than actually needed. 
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III. THE DIFFRACTION PROBLEM 


As noted in the Chapters I and II, one of the basic tasks 
of this thesis research project was to simulate the behavior 
of the proposed multiplexed imaging scheme at long 
wavelengths. Therefore, this first phase of the research is 
called the diffraction problem. It seeks to quantify the 
diffraction behavior of the encoding masks by means of 


computer simulation. 


A. BASIC SCALAR DIFFRACTION THEORY 

No discussion of an imaging system is complete without 
considering the effects of interference and diffraction of the 
incident wavefronts. The importance of these effects cannot 
be understated. 

"Diffraction effects are accordingly oof great 
significance in the detailed understanding of devices 
containing lenses, stops, source slits, mirrors, and so 
on. If all defects in a lens system were removed, the 
ultimate sharpness of the image would be limited by 
diffraction." (Ref. 3] 

The essential feature of diffraction is that an obstruction 
deviates the rectilinear propagation of a light wave. The 
geometry of diffraction through a rectangular aperture is 
Summarized by Hecht [Ref. 3] and depicted in Figure 3. 


References to the aperture plane are in lower case letters; 


those to the image plane are in capitals. 
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Aperture Plane 


Fourier Plane 


Figure 3 Diffraction Through a Rectangular Aperture 

As predicted by the MHuygens-Fresnel principle, a 
differential element of surface area within the aperture, dS 
(equivalent to dydz), may be considered as comprised of an 
infinite number of secondary wave sources. However, for a 
source of monochromatic light, dS is much smaller than the 
wavelength A, and all of dS'‘s contributions to the 
electromagnetic disturbance at point P in the image plane 
arrive in phase and interfere constructively. Assuming 
uniform, normally incident illumination over the entire 


aperture, it follows that 
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E( Ky pi) ee | | 4G) e ~iK(ysing + zsin8) dydz , 


(5) 


where the aperture function is 


a 1 for (y,z) in transparent regions 
(Y,2) =}6 for (y,z) in opaque regions : 
(6) 


The wave number x, is 27/A. The spatial frequencies in the y 


and z directions are, respectively, 


ky = ®Y/R = «sing , 
(7) 


and 
k, = *Z/R=xsin® . 
(8) 


E represents the diffracted electric field distribution. 
.Figure 4 illustrates the angular relationship of @ and @ 
between the aperture and image planes. 

Equation (5) indicates that the electromagnetic field 
distribution at the image plane is proportional to the two- 
dimensional Fourier transform of the aperture function. 
Determining the diffraction behavior of the photon flux 
through the Walsh masks requires a three-step calculation. 
First, simulate the aperture function for each mask of 
interest. Second, calculate the two-dimensional Fourier 
transform of that mask. Third, calculate the actual 


distribution of the resulting irradiance as a function of 
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Figure 4 Diffraction Geometry 


spatial frequencies k, and k,. The third step is simple; the 
irradiance I is proportional to the time average of |E|° 


[Ref. 3), such that 


I « ( EE*) 
(9) 


B. GENERATION OF WALSH MASK APERTURE FUNCTIONS 

The orthogonal Walsh aperture functions were simulated by 
means of ae straightforward algorithm, which is’ easily 
described via an analogy with the familiar topic of Fourier 
transforms. To recover one of the orthogonal trigonometric 


basis functions upon which Fourier analysis is built, one 
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could, in principle, calculate the Fourier transform of a 


Dirac delta function pair [Ref. 5], such that 


| €8 (4-40) + 6(f£+f5) Je 


— 


27Tifpt ~Tifot 2rifot 
0 ere 0 


df-e 
= 2 cos(27fj,t) 


(10) 
While there are more efficient methods of calculating a simple 


cosine function, the transform approach is conceptually sound. 
Similarly, if one has at their disposal an efficient Walsh 
transform routine, it becomes a trivial matter to generate 
Walsh functions of different sequencies: simply Walsh 
transform a numerical equivalent of a delta function at the 
desired sequency. 

The mathematical construction of the Walsh masks as 
outlined above was accomplished by Davis using a fast Walsh 
transform algorithm which he developed from the generalized 
transform algorithm theory presented in Elliot and Rao [Ref. 
7 This program was originally written in Microsoft C. 
Davis translated the code into FORTRAN for incorporation into 
the mask diffraction program WMASK.for (listed in Appendix A) 
aS a subroutine called FWT. The necessary parameters that 
govern the mask generation and ultimate diffraction pattern 
are the desired spatial sequency, the size of the basis set of 
Walsh functions, and the overall dimensions of the two- 
dimensional mask matrix which will be produced. For this 


stage of the project, it was decided that the number of 
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sequencies in the basis set would remain fixed at 64, 
corresponding to an eight-by-eight array of square pixels, as 
shown in Figure 2. The interactive portion of program 
WMASK.for requests specification of the desired mask sequency 
in the range zero to 63. 

The FWT subroutine operates on a delta function and 
returns a one-dimensional array of 64 elements corresponding 
to the discretely sampled Walsh function with the specified 
sequency number. These 64 elements have assigned values of +1 
or -1, representing the transparent and opaque regions, 
respectively, of the Walsh masks. The program adds the value 
one to each returned element, and then divides each element by 
two. This normalizes the sequency array to values of either 
one (for a transparent pixel) or zero (for an opaque pixel). 
As a visual check, whose importance was later realized, the 
program writes the 64 elements to the computer screen in the 
form of a eight-by-eight matrix. Trial and error simulations 
and numerical experiments revealed that this basic eight-by- 
eight matrix format was undersampled for subsequent Fourier 
transformation (Ref. 5], so the matrix structure was expanded 
by a factor of two in each dimension. For clarification, the 
following illustrates the process in generating the Walsh mask 
for sequency 28 of 63 (the 16-square checkerboard in Figure 
2). The 64 element array returned by the subroutine FWT is 


ordered into an eight-by-eight matrix in the form 
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1.1,.0,0,1 Teero 
i l°O0 0 T fey 
00110011 
00110011 
110031100 
1100313100 
00110011 
0 O21 12050 1.2 


After the bi-directional expansion by a factor of two, the 
eight-by-eight matrix becomes a 16 by 16 matrix appearing on 


screen as 


DODORPKFPFRFPOOOOFPFHPPHH 
CODOPFPHRFPOOOOFPEHPFHEH 
DODOPKFPFRPOOOORPKRFRPH 
CODOPRFPHFRPOOOOFRPKREHPH 
PRR rPODOOORPFHFRFHOOOO 
PRR rRPODOOORPPHFHOOOO 
PrHRPrHFODOOORPPFHFOOOO 
PRR HODOOORPFPPFOOOO 
SCODOOPFPHHFPOOOOPFRHPH 
DOCOPFPHHOOOORPFHPH 
DODOOHPKRFHHFPOOOORPFRHPH 
DOOOHPRFPHFFHFPOOOOPFRPP 
PHRHErPOOOORPPRFPHOOOO 
BPHEErPOOOORPHHHOOOO 
PRR HODOOORPPFHPOOOO 
BPE HPHOOOORPHPHEHOOOO 


The two-dimensional Fourier transform subroutine used in 
the program WMASK.for mandates that the input matrix contain 
complex numbers. The required format calls for the real part 
of the number followed by the imaginary part, in column order. 
However, the elements which comprise the Walsh mask are real 
numbers. For this reason, the numerical model of the Walsh 
mask as shown above must be modified by injecting a row of 


zeros, representing the imaginary component of each element, 
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between each row. This transforms the mask to the proper 
complex representation. The resulting matrix is 32 by 16. 

Finally, the mask was centered on a 64 by 32 zero matrix 
in complex form, which constitutes a border of zeros, or 
zeropad, of a factor of two surrounding the mask. This 
dimension doubling was again reached by trial and error. Its 
purpose was to interpolate the resulting spectrum via 
oversampling (Ref. 5]. By increasing the overall dimensions 
of the input matrix by a factor of two (reaching final 
dimensions of 32 by 32), the inherent spectral sampling rate 
is doubled, and consequently, the resulting’ spatial 
frequencies in each dimension of the diffraction pattern are 
then more densely spaced by a factor of two. 

As outlined in Press, et al., [{Ref. 6], the complex data 
must first be properly ordered prior to performing the two- 
dimensional Fourier transform. This ordering is required to 
reduce the number of steps necessary to complete the 
transformation. As is typical of many one-dimensional FFT 
programs, the input data must be reordered in a wrap-around 
fashion by placing the second half of the data in front of the 
first half prior to conducting the transform. This quadrant 
reversal is also discussed in Brigham [Ref. 5]. In two 
dimensions, with the data in matrix form, this wrap-around 
procedure must be performed twice, once in each spatial 
direction. Finally, the matrix is passed to the two- 


dimensional FFT subroutine in normal FORTRAN order. This 
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ordering requires partitioning the matrix into a one- 
dimensional array composed of the sequence of the vertical 
columns of the input matrix. Such an ordering structure is 
well represented in Figure 12.11.1 of Press, et al., [Ref. 6]. 
The FFT subroutine also returns the transformed data in this 
order, which must be sorted following the reverse of the input 
procedure, to recover the two-dimensional FFT of the input 
matrix. 

The following illustrates the steps involved in the data 
manipulation of a simple two-by-two input matrix enroute to 
the FFT subroutine. The numbers represent the matrix 
elements. Each step was originally developed as a separate 
program to verify its function, and then inserted into the 


master program WMASK.for as subroutines. 


a 1 2 Input Data Matrix 
3 
2 1 2 Inject Zeros to Transform to Complex Data 
0 O 
3 4 
0 O 
a7 0 0 0 0 Center Complex Matrix on Zeropad 
0 0 O O 
0. cle 2° a@ 
0 0 dO O 
0 3 4 O 
0 0 O O 
0 0 0 0 
0 0 Oa O 
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a3 0 0 OO 0 Horizontal Wrap-Around 
0 0 O QO 
2 0 Q 1 
Os “ORO 
4 0 0 3 
0 0 O QO 
0 0 O O 
0 0 0 O 
S.. 4 Ome 3 Vertical Wrap-Around 
0 0 0 QO 
0 0 OQ O 
0 0 0 Q 
0 0 Q OQ 
0 0 Q QO 
2am ck 
0 0 90 O 
6. 4 0 0 3 Form 1-D Array to Subroutine 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
2 0 ) 1 
0 0 0 0 


In returning the transformed data, the steps are taken in 


reverse. 


C. THE FAST FOURIER TRANSFORM 

The fast Fourier transform (FFT) has played a significant 
role in reducing the amount of computing time required for 
Fourier analysis. In two or more dimensions, as is required 
by this and other imaging projects, the number of calculations 
becomes so large that, without the FFT, it could not be done 
in a useful amount of time. The key feature which makes the 


FFT work so efficiently is an algorithm to reduce redundant 
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arithmetic operations. Brigham (Ref. 5] outlines the theory 
of the FFT and describes its implementation. 

The two-dimensional FFT used in this project is an 
adaptation of one taken from Press, et al., [Ref. 6], which 
was developed by N. M. Brenner of Lincoln Laboratories. The 
program provided in the reference was written in FORTRAN. 
Consequently, all programming for this part of the thesis was 
done in this language. Although several personal computers 
were used throughout the programming phase, the majority of 
the work was performed on a 16-bit Imperial 80286 machine, 
using the Microsoft FORTRAN 4.1. compiler. To minimize 
computing time, an executable file of the successfully tested 
program was installed on an IBM Model 70 PS/2 80386 personal 


computer with an 80387 coprocessor. 


D. CONVERTING THE TRANSFORMED MATRIX INTO DIFFRACTION 

PATTERNS 

The array returned from the FFT subroutine is reordered 
back into complex matrix form, but does not yet represent the 
diffraction pattern of the input Walsh mask. In reference to 
equation (5), the two-dimensional Fourier transform of the 
aperture function is an expression for the electric field 
distribution E(k,,k,) at the Fourier plane. However, the 
desired quantity is the irradiance I, as a function of the 


spatial frequencies. 
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The algorithm for implementing equation (9) is contained 
within the subroutine INTENSITY of program WMASK.for. MThis 
subroutine converts the returned complex field distribution 
Matrix into a real matrix where each element represents the 
intensity at a point in the Fourier diffraction pattern. In 
doing so, the matrix is reduced from dimensions of 64 by 32, 
to a 32 by 32 square. The matrix then represents the field 
intensity as a function of spatial frequency, which is the 
desired form for the diffraction pattern of the Walsh mask 
aperture function. 

A consequence of performing a two-dimensional fast Fourier 
transform on a square matrix is that, in the final form, the 
first row and the first column contain all zeros. This result 
stems from the vertical and horizontal quadrant reversal 
required in the matrix manipulation prior to transformation. 
As this first row and column have no real meaning, they are 
eliminated from the intensity matrix. Therefore, this last 


operation leaves the matrix in the form of a 31 by 31 square. 


E. PLOTTING THE DIFFRACTION PATTERNS 

It was necessary to locate a software package which could 
produce a three-dimensional perspective plot of the intensity 
matrix as a function of spatial frequency. Such a program is 
SURFER, which is distributed by Golden Software, Inc. This 
program was used to produce all of the diffraction pattern 


plots which appear in Appendix B. 
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The format required by SURFER demands that the data to be 
plotted be accessible in three-dimensional Cartesian (x,y,zZ) 
coordinates. At this stage in the programming, however, the 
intensity was represented as a 31 by 31 square matrix. A 
subroutine called MAT2XYZ was developed to convert the matrix 
into the required format. In this form, each converted point 
was assigned a set of (x,y) coordinates to locate the point in 
the intensity matrix. The (Zz) coordinate indicated the 
intensity at that point. Such a process generates 961 data 
points from the 31 by 31 intensity matrix, and writes them to 
a Gisk file for later use by the plotting routine. 

SURFER employs a two-step process to create three- 
dimensional plots of input data. The first step involves 
forming a regularly spaced grid from irregularly spaced data. 
The data, in the case of this project, are the (x,y,2Z) 
elements which are generated by the master program as 
described above. These data points are imported from the disk 
file into SURFER by specification during the gridding phase. 
In the language of the software, this process is simply called 
"gridding", and is invoked in the menu heading GRID. The 
second step in the plotting process calls for loading the 
gridded data into the subprogram called SURF. SURF is a menu- 
driven, interactive graphics routine which produces’ the 
perspective plots constituting the diffraction patterns. Its 


output is a plot file which contains the SURFER-generated 
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surface plot data. This plot file may be viewed on screen or 
sent directly to a specified plotter. 

As noted, SURFER is a menu-driven program. The following 
parameters were used during the interactive set-up to create 
the diffraction pattern plots which appear in Appendix B. The 
parameters are listed in the order in which they were input, 
and are the result of much trial and error to find the optimal 
settings. Many of the default parameters were satisfactory 
for making the diffraction pattern plots; only the parameters 


which differ from the defaults are listed below. 


At the Main Menu: GRID 

Filespec to Pass: XYZ.GRD 

At Menu: Modify 
Select: Smooth 
Select: Spline (expansion factor 2) 
Select: Begin 


This procedure creates the output grid file OUT.GRD 


At the Main Menu: SURF 
At Menu: Input 
Input grid file: OUT.GRD 
At Menu: View 
Projection: Perspective 
At Menu: LineTyp 
Plot constant (x,y,Z) : x,Y 
At Menu: Base 
Plot base: No 
At Menu: Title 
Plot Title Diffraction Pattern 
For Walsh Mask #28 
Title symbol set: DEFAULT. SYM 
Title position: (0,-1) 
Title angle in degrees: 0 
Title character height: 0.15 
Title color: 1 
Plot orientation legend: No 
Legend position: Automatic 
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At Menu: 
Plot 
Axis 


axes: 

color: 

Axis symbol set: 

Axis plane: 

Tic distance: 

Labeled tic frequency 
Label format: 


Number of decimal digits: 


Label character height: 
Tic label angle: 
Axis title: | 
Title character height: 
Title to label distance: 


Post. 


At Menu: 
Name of plot file: 
Scale factor: 
Page position: 
Plot file format: 
File write mode: 


Number of decimal digits: 


Send to output device: 


Axes 
X,Y (each separately) 
z 


DEFAULT. SYM 
Automatic 
6.28 

2 

Fixed 


2 

0.1 

0 

Spatial Frequency K Sin A/B 


0.13 
0.14 


No changes are required in menu headings 8ize, XYLine, or 


Output 
OUT. PLT 

1 

1,1 
Binary 
Overwrite 
3 

No 


No changes are required in menu heading Equip. 


This procedure generates the plot file OUT.PLT, which is 


called for in the final phase of SURFER. 


the HP7550A plotter. 


At Main Menu: 
Name of plot file: 


Do you wish to shift the entire plot? 
Do you wish to scale the entire plot? 
Scale factor in X direction: 
Scale factor in Y direction: 
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This file is sent to 


PLOT 
OUT. PLT 
No 

Yes 

0.6 

0.6 


F. ANALYSIS OF DIFFRACTION PATTERNS 

Having established a working diffraction program and 
stored the listed plotting parameters into SURFER, the 
remainder of this phase of the thesis research was devoted to 
generating the Walsh masks and to examining their diffraction 
pattern plots to reach conclusions about their diffraction 
behavior. Analysis of the first few diffraction pattern plots 
immediately raised questions about what they represented; the 
results did not seem to conform to expectations. The computer 
code was reviewed to determine where, if any, mistakes were 
present. "What should the diffraction patterns of the Walsh 
masks should look like?" became the question of interest. 
Consequently, the project focus shifted to resolve this 


matter. 


G. APPARENT PROBLEMS WITH THE DIFFRACTION PATTERNS 

The first few diffraction patterns plotted were very 
encouraging. In reference again to the Walsh masks in Figure 
2, sequency zero appears in the upper left hand corner as the 
solid black square. When centered on the zeropad, the black 
Square iS seen aS a square aperture (considering the black 
portion as transparent, the zeropad as opaque). The 
diffraction pattern of a square aperture is well known and 
very identifiable. It appears as Figure 10.25(b) in Hecht 
(Ref. 3]. It was with much relief that the plot generated by 


SURFER of Walsh sequency zero (Figure B-1) was identical to 
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this figure within the scaling limits. Continuing across the 
first row of masks, with sequencies one, two, etc., the 
aiffraction patterns bore the features which would have been 
expected for an array of vertical slits. Despite the apparent 
success, the diffraction pattern produced by sequency 63 (the 
64-square checkerboard in Figure 2) did not conform at all to 
expectations. 

This checkerboard pattern has the finest detail of all of 
the masks in the basis set of Figure 2. Consequently, its 
diffraction pattern would be expected to be spread out much 
more than the others. Also, its symmetry about the main 
diagonal predicts a likewise symmetric diffraction pattern. 
The plot generated from sequency 63 was neither well dispersed 
nor symmetric. Only after installing the routine to display 
the mask on the monitor did the explanation for this become 
apparent. 

1. False Assumptions in the Sequency Ordering 

From the outset, it was assumed that the Walsh masks of 
Figure 2 were sequentially ordered left-to-right in a raster 
fashion, beginning at the upper left hand corner. Certainly, 
the black square represents sequency zero. In this system 
then, the masks of this figure appear "numbered" in sequency 


as 
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In such a system, the first row appears in order, and the 
checkerboard is located in position 63. 

The revelation caused by the routine to display the masks 
indicated an entirely different ordering. By calling for the 
masks in ascending order (from zero), it became evident that 
instead of a raster structure, the subroutine FWT returned the 
data in a_ serpentine fashion. Regardless of this 
understanding, the first row appears identically in each 
scheme. This explains why the square aperture pattern could 
be generated without full comprehension of the sequency 
ordering. However, the mistakenly assumed raster ordering 
places the checkerboard in a different location as opposed to 
where this mask properly appears when located with a 
serpentine ordering. Therefore, rather than generating the 
diffraction pattern of the checkerboard by specifying sequency 
63, the diffraction pattern of the mask located in position 56 
was being produced, which is’ predictably neither well 


dispersed nor symmetric. 
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2. The Proper Ordering of the Walsh Masks 
With the preceding discussion in mind, the proper 
"numbering" of the Walsh masks as generated by the diffraction 


pattern program and appearing in Figure 2 is 


Throughout the remainder of this thesis, the Walsh masks and 
their subsequent diffraction patterns will be referred to 
using this system, whereby the assigned number indicates both 


the sequency of the mask and its location in the above table. 


H. SYMMETRIES IN THE WALSH MASKS 

With a complete understanding of the ordering of the 
masks, it was quickly determined that the program had been 
working all along, and production of the diffraction pattern 
plots resumed in earnest. However, not all of the masks were 
used. Collectively, an examination of the masks in Figure 2 
reveals that masks reflected across the main diagonal are 
identical except for a 90 degree rotation. This symmetry 
allows one to draw conclusions about the complete set of masks 
by examining only the upper (or lower) triangle of the masks 
in this figure. Such a symmetry reduces the number of 


required analyses and plots from 64 down to 36, as the 
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diffraction pattern of a rotated mask is the same as its non- 
rotated counterpart; the only difference being that the 
pattern is also rotated by the same amount as the mask. The 
diffraction patterns which comprise Appendix B, then, are 
taken from the 36 masks which appear on and above the main 


diagonal of Figure 2. 


I. A TYPICAL DIFFRACTION PATTERN PLOT 

Figure 5 shows the diffraction pattern for Walsh mask 28. 
This mask is the 16-square checkerboard. This typical pattern 
illustrates the type of information which is analyzed when 
examining these plots. Angles A and B correspond to angles ¢ 
and ® shown in Figure 4; they relate to, respectively, the 
angular spread due to diffraction in the Z and Y directions of 
the Fourier plane. The axes are labeled in normalized units 
of spatial frequency corresponding to k, and k,. Each 
direction is labeled out to values of + 27. The multiples of 
m are used to keep the dimensions consistent in terms of 
spatial frequency. These units of spatial frequency are the 
reciprocal of the Walsh mask pixel size. 

The irradiance (vertical) axis is neither labeled nor 
scaled. This fact is due to the nature of the Walsh masks. 
All masks (except for the square aperture of sequency zero) 
have equal transparent and opaque areas. Therefore, the 
diffraction patterns of all of the masks will transmit the 


same amount of electromagnetic energy, albeit in different 
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Figure 5 Diffraction Pattern for Sequency 28 
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distributions. Hence, the relative irradiance scaling is the 


same for all masks, and absolute scaling is unnecessary. 


J. CONCLUSIONS CONCERNING THE DIFFRACTION LIMIT 

The quantity which concerns us is the spreading of the 
irradiance in each diffraction pattern. Although all of the 
patterns show a characteristic central peak containing the 
majority of the irradiance, some display a significant amount 
of energy spilling out into the periphery of the pattern. For 
reference purposes, Appendix B constitutes a catalog of all 36 
different patterns generated by the Walsh masks in the basis 
set of 64. They may be examined individually to determine the 
limit in terms of the spatial frequency necessary to retain 
the diffracted irradiance; or, perhaps more expeditiously, 
only the extreme cases need be considered. The result of such 
a process constitutes the answer to the research question at 
hand. 

In reference to the Walsh masks in Figure 2, the finest 
detail occurs in mask 56, the dense checkerboard pattern. 
Certainly its diffraction pattern must be considered as a 
limiting case. Nearly all of the masks in the first row show 
diffraction patterns with wide angular dispersion as well, 
particularly masks 1, 2, and 3. Such patterns define the 
extremes of diffraction spreading for this basis set of Walsh 
masks. It is important to ensure that the limit chosen 


encloses nearly all of the significant irradiance of even the 
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most severe case of diffraction. After careful consideration, 
it was determined that this worst case occurs for Walsh 


sequency three (Figure B-4). 


K. FORMULATING THE DIFFRACTION LIMIT 

Having made the above determination, it became a matter of 
simply reading from the plot the minimum spatial frequency 
required to encompass the majority of the diffracted 
irradiance. It was concluded that this minimum spatial 
frequency was 6 reciprocal length units, as defined above. To 
put this result in a form useful in future instrumental 
design, it may be stated that the following relation must be 


satisfied 


[27 
aa 


sinél <6. 
(11) 


The angle © is used here since, once a mask is "rotated", in 
that its counterpart across the main diagonal is used, the 
angle A becomes equivalent to the angle B. For simplicity, 
the diffraction angle henceforth shall be called @. Figure 6 
illustrates the simplified relationship between the mask, the 


optics, and the detector. 
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Mask Optics Detector 





Figure 6 Mask, Optics, and Detector Orientation 
The formulation of equation (11) constitutes the solution 


to the first thesis problen. 
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IV. RESOLUTION CRITERIA FOR IMAGE PATTERN RECOGNITION 


A. OVERVIEW 

The second phase of this research project involved, by 
means of computer simulation, a determination of what minimal 
Spatial resolutions are needed for pattern recognition in 
images of commonplace objects. Establishment of these 
criteria are of great importance in the overall multiplexed 
imaging research project, for they dictate the minimal number 
of picture elements (pixels) needed by a practical remote 
sensing instrument. This minimal number of pixels, in turn, 
specifies the number of Walsh basis functions and encoding 
Masks required to accomplish this task. 

A working prototype system does not yet exist; therefore 
this problem was addressed in reverse. Digitized images were 
transformed into a sets of sequency spectra. These spectra 
were sequency-truncated (low-pass filtered) to mimic the 
effects of limited resolution in an actual instrument. The 
filtered sequency spectra were then re-transformed into the 
image domain for analysis. 

The procedure may be further clarified by recalling that 
an image's higher spatial frequencies are responsible for 
finer image detail. In this research, spatial sequencies play 


the role normally reserved for spatial frequencies with regard 
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to traditional image analysis approaches. Elimination of 
higher sequencies degrades image resolution in the same sense 
that the fine detail in an observed field is degraded by the 


loss of components of higher-order spatial frequency. 


B. THE WALSH FUNCTIONS AND IMAGE PROCESSING 

Previous sections have discussed the Walsh functions with 
regard to sequency theory and the proposed encoding technique. 
Beauchamp [Ref. 4] goes much further in detailing several 
other applications of the Walsh functions in image processing, 
including pattern recognition, and image enhancement, 
transmission, and restoration. However, those properties of 
the Walsh functions which are summarized in Chapter II of this 
thesis are sufficient for a description of the procedures used 


in this phase of the research. 


C. CREATING DIGITIZED IMAGES 
The actual order of tasks for determining the number of 
sequencies required for addressing the image resolution - 


pattern recognition problem were as follows: 


e Capture a series of images on a videocamera 

e Digitize the images into computer disk files 

e Write software to store, access, and display the images 
e Write software to perform the sequency low-pass filtering 


e Analyze the resulting images to establish appropriate 
resolution criteria. 
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D. CAPTURING THE IMAGES 

The original plan was to record still images and to 
perform a digitizing scan of each using a LYNXX solid-state 
CCD slow-scanning imaging camera. Prof. D.L. Walters procured 
such a camera system along with the necessary imaging 
software, from SpectraSource Instruments. However, the camera 
system was modified for higher frame rates, and was shipped 
back to the manufacturer. Hence, because the system was not 
available for use in this project, an alternative was sought. 

The most readily available device for image recording is 
an ordinary VHS camcorder. A Panasonic VHS-C Palmcorder was 
borrowed from a friend to record the images used in this 
project. A series of 15 images were recorded from areas of 
interest around the Monterey Peninsula, with the understanding 
that the VHS tape would be played back ina freeze-frame mode 
into digitizing equipment operated by the NPS Physics 
Department. The output from this equipment would then be a 


digitized disk file of each image in simple binary format. 


E. DIGITIZING THE IMAGES FROM VHS TAPE 

Difficulties arose with regard to gaining access to the 
original image digitization equipment, due to the author's 
lack of familiarization with it. Therefore, a quest to locate 
alternate means to complete the task was undertaken. An ideal 
system would grab a frame from the VHS playback, scan it, 


digitize it, and then write the image to an binary disk file. 
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Several systems were located throughout the NPS campus which 
apparently had this capability. 

1. Apple Macintosh IIx Programs 

Two programs in use at the NPS Linear Accelerator facility 
were offered by Prof. X.K. Maruyama, and each was thoroughly 
investigated for suitability to the image digitization task. 
The programs are IMAGE 1.29, and PIXELWORKS. Both programs 
run on an Apple Macintosh IIx personal computer, and have 
Similar capabilities well beyond the simple needs of this 
project. The programs are able to manipulate and digitize an 
image, which is played back from the camcorder into the video- 
in port of the Macintosh. Difficulties arose at this stage of 
the process. 

The programs offer a wide range of digitized data formats; 
in particular, TIFF (and simple TIFF), PICT, MacPaint, 
Palette, and RAW, are available. However; binary is not. 
Substantial effort was made to resolve this problem, but 
despite a collaborative effort, no simple solution was found. 
As an alternative approach, Davis suggested that if the format 
of the data type was known, he could then write a program to 
convert the data into binary bytes. Unfortunately, the 
supporting documentation for the Macintosh imaging software 
was not detailed enough to provide this information. 
Therefore, further investigation into this translation method 


were halted, and alternate digitization means were sought. 
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2. Other Means of Digitization 

Earlier failures in producing a digitized image disk file 
in binary format lead to an expanded search for an alternate 
means to accomplish this task. Several Hewlett-Packard page 
scanners are located in different departments at NPS, and were 
considered as possible candidates. The NPS Computer Center 
also uses HP scanners, as well as an IBM page scanning system. 
The IBM system appeared most promising, but it proved to 
demand a critical investment in terms of training time and 
familiarization. Furthermore, it was perpetually 
oversubscribed. Eventually, these page scanners were dropped 
as alternate means of digitization. 

3. Successful Digitization 

Davis contacted Prof. A.W. Cooper and Mr. W.J. Lentz of 
the NPS Physics Department, who agreed to loan some of their 
equipment and expertise to assist the author with the image 
digitization. Their equipment configuration mandated that the 
original VHS<-C camcorder images be translated onto a standard 
VHS cassette tape. This was easily done using the playback 
feature of the camcorder in conjunction with the record mode 
of the author's VHS recorder. 

The digitizing equipment made use of integrated software 
coupled through an IBM personal computer to link a playback 
VCR and a display monitor. The digitizing program was DT - 


IRIStutor v01.01, which is marketed by DT, Inc. This program 
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utilizes an IRIStutor DT2861 frame-grabber board . in 
conjunction with the IBM PC to reduce a captured image to a 
512 by 512 pixel digitized image file, in eight bit binary 
format. The images from the VHS tape were subsequently 
Gigitized and stored on 5-1/4 inch 360k byte diskettes, with 


one image per diskette. 


F. CONVERTING AND DISPLAYING THE IMAGES ON SCREEN 

The files generated by DT-IRIS contained a header block of 
512 bytes prior to the actual pixel-by-pixel digitized 
representation of the image. It was necessary to strip off 
this information prior to displaying the image on the VGA 
monitor of the project's IBM PS/2 model 70 PC. In order to 
accomplish this, Davis wrote a program in Microsoft C5.1 
called C512, which is used to convert the 512 by 512 DT-IRIS 
images into simple eight bit binary format. This program is 
included as Appendix J. 

Once the image file had been converted into usable forn, 
it was necessary to write another routine to display the image 
on the monitor. Again programming in C, Davis developed the 
program D512, which accesses the binary image file, and 
displays it in 16-level greyscale on the VGA monitor. In this 
scale, a shade of zero represents black, while a shade of 15 
is bright white. To accommodate the line restrictions imposed 
by the monitor, the program reduces the display to 480 rows by 


512 columns. This is accomplished by stripping off the top 
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16, and the bottom 16, horizontal lines of the image. Thus, 
in the absence of any segquency filtering, the program displays 
the image in 480 by 512 pixel format. The program D512 is 


included as Appendix K. 


G. FIVE REPRESENTATIVE IMAGES 

Originally, a total of 15 scenes were recorded using the 
Panasonic camcorder. These 15 sane were viewed with Davis 
to select those which contained characteristic features and/or 
desirable symmetries which might prove interesting in the 
context of the low-pass sequency filtering problen. The 


following images were chosen: 


e U.S. Coast Guard Cutter Point Barrow 

e U.S. Coast Guard Lighthouse at Point Pinos 
e A friendly German shepherd 

e A house with a picket fence 


e U.S. Coast Guard emblem 


Despite the apparent Coast Guard bias in subject material, 
the images selected represent a broad range of spatial 
characteristics such as fine detail, periodic structure, sharp 
contrast, circular symmetry, and so forth. The original, 
unprocessed, images appear in the Figures 7 through 11. Note 
that each figure is captioned with a noun descriptor followed 


by a number. This naming convention has been adopted 
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throughout the remainder of the research. The noun name 
simply describes the dominant object in the image; the number 
indicates the number of sequencies used in each direction in 
the image basis set. Using this scheme, the first image is 
Cutter512, indicating that it is the Coast Guard cutter with 
no sequencies eliminated (i.e., 512 by 512 pixels). 

This labeling system is used throughout the project and is 
most useful when studying the images contained within 
Appendices C - I. However, there is one departure from the 
systen. Images labeled with a "V" or an "H" postscript 
indicate that the image was generated using only a one- 
dimensional Walsh transform, with that transform applied in 
the vertical (V) or horizontal (H) direction. These 
descriptors appear only in conjunction with the sequency 
filtering applied to the image of the cutter, and are 
contained in Appendices C and D. 

All of the images which appear in the appendices were 
photographed directly off of the VGA monitor. For this 
reason, the curvature of the screen is noticeable. The small 
black tab which seems to move about on the left and lower 
borders is a marker placed by the author to identify the 
developed prints. The shaded diagonal strip which is evident 
in some of the photographs, such as in Cutter2V in Figure C-9, 
is a consequence of the monitor's 1/30 second refresh and 
flyback time when photographed. The pictures were taken with 


black and white ISO 400 KODAK Tri-X pan film, at an aperture 
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setting of £4, and shutter speed of 1/15 second. The camera 


was a Canon AE-1. 


H. PROCEDURES BY WHICH IMAGE RESOLUTION WAS MODIFIED 

The FORTRAN program IMAGE.for, written by the author, is 
the primary program which Simulates image resolution 
degradation. It does this by taking an original, full 
resolution, digitized image produced by program C512, and 
operates on it to generate a new file with degraded 
resolution. The degradation involves three steps. First, the 
program calculates the Walsh transform of the full-resolution 
image, to produce its sequency spectrum. Second, it performs 
low-pass sequency filtering on the image spectrum, with the 
maximum sequency cutoff specified by the user. Third, it 
calculates the inverse Walsh transform of the filtered 
spectrum, yielding a modified image. Each of these steps is 
described in detail below. The reader may also wish to 
consult Appendix L, which contains a listing of program 
IMAGE.for. 

1. Production of the Image's Sequency Spectrum 

The program IMAGE.for uses the same fast Walsh transform 
Subroutine (see subroutine FWT in Appendix L) as_ was 
incorporated into program WMASK.for, described earlier in 
section III.B. of this thesis. Two different types of 
transform geometries were used in this research. For most of 


the images of interest, a straightforward two-dimensional 
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Figure 7 Cutter512 
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Figure 8 Lhouse512 
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Figure 9 Dog512 
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Figure 10 Fence512 
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Figure 11 Emblem512 


on 


resolution degradation was required, with the same sequency 
filter cutoff in both X and Y directions. A second approach, 
which was applied only to the Cutter image, investigated the 
effects of one-dimensional sequency filtering along both the 
X and Y directions, independently, while leaving the other 
directions unfiltered. 

The one-dimensional transform procedure is the simpler of 
the two. If resolution degradation is desired in the xX 
(horizontal) direction, the image is first partioned into 
rows. The program calls the FWT subroutine for that row, and 
returns its sequency spectrum. It then cycles to the next 
row, and loops until all rows' spectra have been computed. 
The one-dimensional transform for resolution degradation in 
the Y (vertical) direction is effectively the same, except 
that the program operates on columns, rather than rows. 

The two-dimensional transform operation is essentially a 
combination of both one-dimensional versions. The program 
first performs the Walsh transform on each column of the image 
in turn, then it calculates the transform for each row. This 
ylelds the two-dimensional sequency spectrum of the original 
spectrum. Clearly, the efficient implementation of such a 
numerically intensive calculation is completely dependent upon 
a fast transform algorithm; direct calculation would be 
prohibitively slow. Even with the FWT, full two-dimensional 


Walsh spectrum calculations for a 512 by 512 pixel image 
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required nearly 300 seconds on a 25 MHz2 80386 PC, using an 
80387 math coprocessor. 

2. The Sequency Filtering Operation 

When the program IMAGE.for is first invoked, it prompts 
the user for sequency filter cutoff information. The actual 
sequency filtering operation is then quite simple. The 
program cycles through each sequency spectrum row and/or 
column that is to be filtered, and sets spectral samples at 
sequencies beyond the cutoff sequency to zero. This is not 
unlike the zeropad procedure described in section III.B., 
although in this case the size of the spectrum array is not 
increased. 

3. Calculation of the Degraded Image 

The inversion of the filtered sequency spectrum to its 
corresponding image is straightforward: IMAGE.for calls FWT to 
calculate the inverse Walsh transform of each row and/or 
column of the spectrum to yield the image. Fortunately, this 
is a particularly simple procedure using the type of algorithm 
developed by Davis for subroutine FWT. When the data are 
retained in increasing sequency order (as they are here), and 
the binary bit-reversal sorting portion of the fast transform 
is followed by a decimation in sequency, the Walsh transform 
and its inverse are algorithmically identical, except for a 


simple multiplicative scaling factor following the inverse 


5) 


transform operation. This property greatly simplified the 


writing of program IMAGE.for 


I. ANALYSIS OF RESULTS OF THE RESOLUTION DEGRADATION PROCESS 

A complete set of resolution degradation results for each 
of the five original images is included in its own appendix. 
The reason for this is that the conclusions which are drawn 
from these series of image photographs are highly subjective; 
it is not expected that all readers of this thesis will agree 
completely with the interpretations. Each series of sequency 
filtered degradation is represented in it's entirety for the 
reader to confirm or refute the author's analysis. For 
l1llustrative purposes, the image of the Coast Guard cutter was 
degraded in three separate manners; once vertically, once 
horizontally, and once in two dimensions. All other images 
were degraded in two dimensions only. The appendices contain 
the following series of images and should be consulted 


throughout the following discussion. 


e Appendix C - Cutter(512 O)V 


e Appendix D - Cutter(512 - 0)H 


e Appendix E - Cutter(512 - 0) 
e Appendix F - Lhouse(512 - 0) 
e Appendix G - Dog(512 - 0) 

e Appendix H - Fence(512 - 0) 


e Appendix I - Emblem(512 - 0) 
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1. One-Dimensional Sequency Filtering 
a. Vertical Filtering of the Cutter Image 

Referring to the images of the cutter in Appendix 
C, the effects of vertical sequency filtering are clearly 
illustrated. Figure C-1 is the unfiltered image of the 
cutter, accordingly titled Cutter512. In each successive 
figure, the vertical image resolution is made coarser by a 
factor of two. For example, Figure C-2 represents Cutter256V, 
and is produced from half as many sequencies as Cutter512. 
Figure C-3 represents Cutter128V, which is produced from half 
aS many sequencies as Cutter256V, and so forth. This pattern 
of resolution reduction (in both one and two dimensions) 
continues throughout this and all subsequent series of image 
degradations. 

The basic structure of Cutter512 in Figure C-1 
indicates four vertically spaced regions corresponding to the 
sky, the hillside, the cutter and its moorings, and the water 
in the foreground. With such an orientation of dominant 
features, one would expect that vertical filtering would 
preserve these regions until only the lowest vertical sequency 
components remained. This, in fact, is the case. In Figures 
C-1 through C-5, corresponding to vertical resolutions of 512, 
256, 128, 64, and 32 pixels, this four-region structure 


remains obvious. It begins to become unidentifiable only at 
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vertical resolutions of 16 or less, whereas the cutter itself 
loses recognizable form at a sequency cutoff of 64. 

A noticeable feature of this and the following image 
degradation series is the apparent white "spots" which seem to 
outline the edges of the image. Cutter256 in Figure C-2 
exhibits this behavior quite clearly. It is hypothesized that 
these spots are the result of the type of "ringing", or 
overshoot phenomenon, associated with the processing of 
adjacent areas of high contrast. If this is the correct 
explanation for their presence, then there are more than 
likely just as many black "spots"; the eye just does not see 
them standing out as clearly. The nature of these spots and 
their effect on the imaging process are offered as recommended 
areas of future research. 

b. Horizontal Filtering of the Cutter Image 

Appendix D contains the horizontally degraded 
series of images of the cutter. As before, the first image in 
the series, Cutter512 in Figure D-1, represents the image in 
its undegraded form. All subsequent images in this series are 
degraded by a factor of two in the horizontal direction. 

Cutter512 displays much greater spatial variation 
in the horizontal as compared to the four dominant vertical 
regions discussed above. Therefore, one would expect that 
eliminating the higher order sequencies in the horizontal 


would have less effect on the overall image degradation. This 
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expectation is validated by comparing Cutter256H in Figure D-2 
with the original, undegraded image. When making this same 
comparison with Cutter256V (Figure C-2), it becomes readily 
apparent that horizontal filtering detracts less from the 
image quality (as opposed to vertical filtering) for this 
particular image, due to it's spatial characteristics. 

Another easily seen result of horizontal filtering 
of the cutter image is the rapid degradation of the narrow, 
vertical features. The cutter's mast, for example, may be 
described as such, and is produced from higher order 
sequencies. One would then expect that the mast would quickly 
lose its detail while the image underwent horizontal sequency 
filtering. This is indeed the case. The mast becomes largely 
indistinguishable in Cutterl128H (Figure D-3), while the 
overall image of the ship is still recognizable in Cutter32H 
(Figure D-5). 

2. Two-Dimensional Sequency Filtering 
a. The Coast Guard Cutter 

Here, and for all subsequent image. series, 
resolution is diminished by a factor of two in both vertical 
and horizontal directions for each successive degradation. The 
two-dimensional filtering of this image is displayed in 
Appendix E, and is essentially the combination of the one- 
dimensional horizontal filtering as described above and the 


one-dimensional vertical filtering. Analysis of this series 
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concludes that although the four distinct vertical regions 
(the sky, hillside, moorings, and water) are still visible 
beyond a sequency cutoff of 64, the spatial structure of the 
cutter itself in Figure E-4 is no longer recognizable. 
b. The Lighthouse 

The series of degradations of this image appear in 
Appendix F. The undegraded image, Lhouse512 in Figure F-1, 
seems partitioned into two dominant vertical regions, and 
three horizontal. However, none of these regions incorporate 
any fine detail in either dimension. Therefore, one would 
expect that the image pattern would remain recognizable even 
with comparatively few low-order sequencies. Lhouse64 in 
Figure F-4 confirms this, as the lighthouse is clearly 
recognizable when formed from 64 resolution elements. In 
fact, Lhouse32 in Figure F-5, is still largely identifiable 
from only half again as many elements. 

c. The German Shepherd 

This series of degradations appears in Appendix G. 
The undegraded image in Figure G-1 indicates that there are 
few, if any, right angles in the scene. This figure also 
shows a great amount of fine detail, as well as a somewhat 
circular symmetry around the dog's head. Therefore, one would 
expect that the pattern would require most of the higher order 
sequencies to remain identifiable. This is not so ! Facial 


characteristics, such as eyes, nose, and mouth, appear 
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discernible even with resolution degraded to 32 sequencies as 
seen in Dog32 in Figure G-5. There is most likely a strong 
psychological reason for this. Perhaps the visual information 
processing areas of the human brain selectively register this 
pattern as "wolf!" or "danger!", even with very limited 
spatial image information. 
d. The House with the Picket Fence 

This image was chosen for analysis because of the 
pronounced periodic structure exhibited by both the fence and 
the support columns on the front porch. This is evident from 
Figure H-1. Both features are vertical structures of high 
contrast, and should be susceptible to modification by 
horizontal sequency filtering. The spatial period of the 
columns is greater than that of the fence. Therefore, one 
would predict that the fence detail would be filtered out 
first. Fence256 and Fencei28, in Figures H-2 and H-3, 
respectively, bear this prediction out. Nonetheless, the 
object of primary concern with regard to pattern recognition 
is the house itself, and it loses recognizable form beyond 
resolution with 64 elements. 

e. The Coast Guard Emblem 

This image was selected for analysis because it 
represents several types of symmetry whose degradation effects 
are difficult to predict. The undegraded image in Figure I-1 


shows some horizontally varying structure in the bars of the 
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shield, but the overall symmetry is circular. Otherwise, 
there is very little purely horizontal or purely vertical 
detail, so that one is left to empirical judgment to evaluate 
results. 

Emblem64 in Figure I-4 contributes to the 
conclusion that the shield structure becomes unrecognizable 
beyond a sequency cutoff of 64. However, the bulk circular 
outline of the anchors and lifering are still apparent; in 
fact, they remain quite visible in Emblem32 in Figure I-5. 
This is somewhat surprising because these high contrast, high 
sequency features do not suffer rapid image degradation while 
undergoing sequency filtering. As was the case with the image 
of the dog, results did not agree closely with predictions. 
The common feature between these two departures from 
prediction is the circular symmetry. Recall that the Walsh 
functions which were used throughout this project were 
generated in rectangular coordinates. Perhaps there is a 
connection to the coordinate system of the applied Walsh 
transform and the dominant symmetry of the observed scene. 


Further research into this apparent consequence is in order. 


J. DETERMINING THE PATTERN RECOGNITION RESOLUTION LIMIT 
From the sequency filtered degradations of the five 

images, the author concludes that, with the exception of the 

Emblem series, the images require at least 64 sequencies for 


minimum pattern recognition. Again, it must emphasized that 


60 


this is a somewhat subjective judgment. The Emblem series is 
excluded from this general statement because, as the reader 
should concur, it remains distinguishable beyond this sequency 
cutoff, down to a minimum of 32. It is believed that this is 
related to the highly circular nature of this particular 
image. 

Nonetheless, it is with confidence that this phase 
of the research concludes that the minimum number of 
sequencies required for threshold pattern recognition is 64. 
This result allows one to quantify several physical 
parameters. For example, this sequency limit mandates that 64 
sequencies be available to form the basis set of orthogonal 
functions in each of the X and Y spatial directions. 
Furthermore, the number of pixels required to display the 
output image must be fixed at 4096, which represents the 


Square of 64. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The intent of the initial task of this project was to 
mathematically model the two-dimensional sequency-ordered 
Walsh functions, and to determine the diffraction pattern of 
each. This was accomplished. The result of studying the 
patterns lead to the formulation of the following equation 
which quantifies the diffraction limited performance of the 


encoding scheme: 


aus 
= 


sinél<6. 


The angle © represents the maximum tolerable diffraction angle 
which will ensure that all of the diffracted photon flux 
passing through the Walsh mask will be focused onto the 
detector. This angle, coupled with the size of the detector 
and the wavelength of operation, provides useful information 
to the follow-on work which will proceed to design the optics 
for a prototype instrument. This simple relationship is 
illustrated in Figure 6. 

The second question addressed by this thesis was directed 
towards establishing the minimum number of sequencies required 
to produce a recognizable image. This was accomplished by 
generating images with increasing numbers of higher order 


sequencies filtered out. Once the image pattern could no 
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longer be deemed identifiable, it was concluded that the 
minimum pattern recognition limit had been reached. From 
examining the all of the series of sequency filtered images in 
Appendices E - I, it was concluded that no less than 64 
sequencies are required to produce a recognizable pattern. 

It is recommended that future thesis work be directed 
toward investigating how the coordinate system of the Walsh 
functions affect the results. This question is raised from 
analysis of the series of degraded images of the circular 
Coast Guard emblem (Appendix I). It appears that perhaps the 
dominant symmetry of the object plays an important role in 
determining the minimum number of sequencies required to 
resolve it. An interesting project may be to regenerate the 
findings of this project while using Walsh functions derived 
under a different coordinate system, and make comparisons 
between the two. 

Furthermore, queries into the dependency of image 
resolution as a function of the available greyscale level, and 
the cause and effects of the "ringing" phenomenon apparently 
caused by processing adjacent elements of high contrast, are 
recommended. Either aspect could possibly affect the findings 
of this thesis, and may mislead future research based upon the 


conclusions drawn herein. 
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APPENDIX A 


LT B .J. MUSSELMAN 


THESIS PROGRAM - WMASK.FOR 


THIS PROGRAM COMPUTES MATHEMATICAL MODELS OF WALSH SEQUENCY 
MASKS, AND CALCULATES THEIR DIFFRACTION PATTERNS THROUGH 

2 - D FOURIER TRANSFORMS. THE COMPLEX TRANSFORM MATRIX IS 
REDUCED TO AN INTENSITY MATRIX, WHICH IS STORED IN CARTESIAN 
COORDINATES (X,Y,Z) IN A FILE CALLED C:\FORTRAN\ TEMP\ XYZ. DAT 


THE INPUT MATRIX CWALSH MASK) IS 16 X 16, AND IS CENTERED 
IN A 32 X 32 ZEROPAD. THE ENTIRE MATRIX IS SENT TO THE 
2 - D FOURIER TRANSFORM SUBROUTINE. 


INTEGER ISIGN, NDIM, ROW, NROW, COL, NCOL, NN,N,R, J, SEQ, NUMSAMPLS 
DIMENSION NN (2) 

REAL IDATA(64, 32), TDATA(2048) , PDATA( 1024) 

REAL W(64),X(961), Y(961),Z(961), TEMP, EESTAR (32, 32) 

REAL A(8,8),B(8,16),C(16, 16),CLIN(256) 

NDIM = 2 


NN (1) REPRESENTS ROWS, NN (2) REPRESENTS COLUMNS 


NNC1) = 32 
NNC2) = 32 


32 
32 


NCOL 
NROW 


ISIGN = 1 
NUMSAMPLS = 64 


OPEN (11,FILE = ’°Cs\FORTRAN\ TEMP\XYZ.DAT’, STATUS = ’OLD’ ) 


WRITE C#,#) "WMASK.FOR IS CURRENTLY CONFIGURED TO PRODUCE A BASIS’ 
WRITEC#,#) "SET OF 64 WALSH SEQUENCIES NUMBERED FROM O TO 63.’ 
WRITE C#, #) 

WRITEC#,#) "ENTER SEQUENCY NUMBER TO GENERATE AND TRANSFORM :’ 
READ(C#,#) SEQ 


DO S2 ROW = 1, NROW * 2 
DO Sz COL = 1, NCOL 
IDATACROW,COL) = 0.0 


CONT INUE 

DO 10 J = 1, NUMSAMPLS 
WCJ) = 0.0 

WCSEQO + 1) = 1.0 


CALL FWT CNUMSAMPLS,W,.FALSE. ) 


THIS SECTION DEVELOPS THE 2X EXPANDED MATRIX USING WALSH 
SEQUENCIES AS RETURNED FROM SUBROUTINE FWT. THE MATRIX 
IS WRITTEN TO THE SCREEN AS A VISUAL CHECK. 


J = 1 


DO 32 ROW = 1, SQRT CNUMSAMPLS) 
DO 32 COL = 1, SQRTCNUMSAMPLS) 


64 


“ex & & 


rs 


32 


43 


34 


35 


36 


100 
101 
102 
201 


ACROW,COL) = .5 * (W(J) + 1) 
J=eJ+i1 
CONTINUE 
DO 43 ROW = 1, SORT (NUMSAMPLS) 
DO 43 COL = 1, SQRT(NUMSAMPLS) 
BCROW, 2 * COL) = ACROW,COL) 
BCROW, 2 * COL - 1) = BCROW, 2 #* COL) 
CONT INUE 
DO 34 ROW = 1, SQRTCNUMSAMPLS) 
DO 34 COL = 1, 2 * SQRT(NUMSAMPLS) 
C(2 * ROW , COL) = BCROW,COL) 
C(2 * ROW - 1,COL) = C(2 * ROW,COL) 
CONTINUE 
DO 35 ROW = 1, 2 * SQRT(NUMSAMPLS) 
WRITEC#,201) (CCROW,COL) , COL = 1, 2 * SQRT(NUMSAMPLS)) 
CONT INUE 
READ(#, *) 


Sat 
DO 36 ROW = 1, 2 * SQRT(NUMSAMPLS) 

DO 36 COL = 1, 2 * SQRT(NUMSAMPLS) 

CLINCJ) = CC(ROW, COL) 
J=eJ+1 

CONT INUE 
FORM THE COMPLEX MATRIX WITH ALL IMAGINARY COMPONENTS ZERO 
CALL ZEROPAD (CLIN, IDATA) 
CALL FLIPFLOP(IDATA, NROW, NCOL) 
CALL MAT2LIN (NROW,NCOL, IDATA, TDATA) 
CALL FOURN (TDATA,NN,NDIM, ISIGN) 
CALL LIN2MAT C(NROW,NCOL, IDATA, TDATA) 
CALL FLIPFLOP CIDATA, NROW,NCOL) 
CALL INTENSITY (NROW, NCOL, IDATA, EESTAR) 
CALL MAT2XYZ (NROW,NCOL, EESTAR, X,Y, Z) 
WRITE XYZ MATRIX TO FILE FOR PLOTTING 
DO 15 N= 1, (NROW - 1) * CNCOL - 1) 

WRITE(11,102) X(N), YON), ZCOND 
CONT INUE 
FORMAT(1X,FS. 2) 
FORMAT (1X,8(F8. 2, 2X)) 
FORMAT (1X,F8.3, 2X,FS8. 3, 2X, E12. 6) 
FORMAT (1X, 16(F4.1,1X)) 
END 
SUBROUTINE ZEROPAD (CLIN, IDATA) 


INTEGER J,ROW, COL 
REAL CLIN(1024), IDATAC64, 32) 
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J=1 
DO 12 ROW = 17,47,2 
DO 12 COL = 9,24 

IDATACROW,COL) = CLIN(J) 
JeJ+i 

CONT INUE 

RETURN 

END 


SUBROUTINE MAT2LIN (CNROW,NCOL, IDATA, TDATA) 


INTEGER N,ROW, COL 
REAL IDATAC64, 32), TDATAC 2048) 


N= 1 
DO 4 COL = 1, NCOL 
DO 4 ROW = 1, NROW #* 2 
TDATACN) = IDATACROW, COL? 
N=N +1 
CONT INUE 
RETURN 
END 


SUBROUTINE LIN2MAT (CNROW,NCOL, IDATA, TDATA) 


INTEGER N,COL, ROW 
REAL IDATAC64, 32), TDATAC2048) 


N= 1 

DO 6 COL = 1, NOL 
DO 6 ROW = 1, NROW * 2 
IDATACROW,COL) = TDATACN) 
N=N +1 

CONT INUE 

RETURN 

END 


SUBROUTINE INTENSITY (NROW,NCOL, IDATA, EESTAR) 


INTEGER N,ROW, COL 
REAL IDATAC64, 32) ,PDATAC 1024) ,EESTAR(32, 32) 


N= 1 
DO 68 COL = 1, NCOL 
DO 68 ROW = 1, (2 * NROW) -1, 2 
PDATACN) = CIDATACROW,COL))##2 + CIDATACROW + 1,COL)II##2 
N=N + 1 
CONT INUE 


N= 1 
DO 9 COL = 1, NCOL 
DO 9 ROW = 1, NROW 

EESTAR(ROW,COL) = PDATACN) 
N=N +1 

CONTINUE 

RETURN 

END 


SUBROUTINE MAT2XYZ CNROW,NCOL, EESTAR, X,Y, Z) 
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INTEGER N, ROW, COL 
REAL EESTAR(32,32),X(961),Y(961),2(961) 


N= 1 
DO 2 ROW = 1, NROW - 1 
DO z COL = 1, NCOL - 1 

X(N) = -6.2831 + (CROW - 1) * (€12.5662/30) 
YCN) = -6.2631 + (COL - 1) * (€12.5662/30) 
Z(N) = EESTARCROW + 1,COL + 1) 
N=N + 1 

CONTINUE 

RETURN 

END 


SUBROUTINE FLIPFLOP (BOX,NROW, NCOL) 


INTEGER F,C,NROW,NCOL 
REAL BOX (64,32) 


RE-ORDER MATRIX LEFT TO RIGHT 


DO 2 C = 1, NCOL/2 
DO 2 R = 1, NROW * 2 
TEMP = BOX(R,C) 
BOX(R,C) = BOX(R,NCOL/2 + C) 
BOX(R,NCOL/2 + C) = TEMP 
CONTINUE 


RE - ORDER MATRIX TOP TO BOTTOM 


DO 3 C = 1, NCOL 
DO 3 R = 1, (NROW#2)/2 
TEMP = ROX(R,C) 
BOX(R,C) = BOX((NROW#2)/2 + R,C) 
BOX ( (NROW#2)/2 + R,C) = TEMP 
CONTINUE 
RETURN 
END 


SUBROUTINE FOURN (DATA,NN,NDIM, ISIGN) 


REAL WR,WI,WPR,WPI,WTEMP, THETA 
DIMENSION NNC(NDIM), DATAC#) 
NTOT = 1 


DO 11 IDIM = 1, NDIM 
NTOT = NTOT # NNCIDIM) 
CONT INUE 
NPREV = 1 
DO 168 IDIM = 1, NDIM 
N = NNCIDIM) 
NREM = NTOT / (N #* NPREV) 


IP1 = 2 * NPREV 
IPZ = IF1 # N 
IP3 = IP2 #* NREM 
IZREV = 1 


DO 14 I2 = 1,1P2,IP1 
IF (12 .LT. I2REV) THEN 


DO 13 Ii T2, I2 + IP1 - 2, 2 
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DO 12 13 #11, IFP3, IF2Z 
ISREV = I2REV + 13 - I12 
TEMPR = DATA (13) 

TEMPI = DATA C(I3 + 1) 
DATA (13) = DATA CISREV) 


DATA (13 + 1) = DATA CISREV + 1) 


DATA CISGREV) = TEMPR 
DATA CISREV + 1) = TEMPI 


CONT INUE 
CONTINUE 
ENDIF 
IBIT = IP2/2 
IF CCIBIT .GE. IP1) .AND. CI2Z2REV .GT. IBIT)) THEN 
IZREV = I2REV - IBIT 
IBIT = IBIT/2 
GO TO 1 
ENDIF 
IZREV = I2REV + IBIT 
CONT INUE 
IFP1 = IP1 
IF CIFP1 .LT. IP2) THEN 
IFP2 = 2 # IFPI 
THETA = ISIGN # 6.2831853071717959D0/ CIFP2/IP1) 
WPR = -2.D0O # DSIN(O.SDO # THETA) ## 2 
WPI = DSIN (THETA) 
WR = 1.D0 
WI = 0.DO 
DO 17 13 =1, IFPI, IP1 
DO 16 11 #® 13, 13 +IPi - 2, 2 
DO 15 12 = 11, IPS, IFP2 
Ki = I2 
K2 = Ki + IFP1 
TEMPR = SNGL (WR) #*DATACK2) —-SNGL CWI) #*DATACK2+1) 
TEMPI = SNGLCWR) #£DATACK241)+SNGL (CWI) #DATACK2) 
DATACK2) = DATACK1) - TEMPR 
DATACK2 + 1) = DATACK1 + 1) - TEMPI 
DATACK1) = DATACK1) + TEMPR 
DATACK1 + 1) = DATACK1 + 1) + TEMPI 
CONT INUE 
CONT INUE 
WTEMP = WR 
WR = WR # WPR - WI # WPI + WR 
WI = WI # WPR + WTEMP # WPI + WI 
CONT INUE 
IFeie Sele Pe 
GO TO 2 
ENDIF 
NPREV = N # NPREV 
CONT INUE 
RETURN 
END 


SUBROUTINE FWT CNUMSAMPLS, DATARRAY, ENABLESCALE) 
FAST WALSH TRANSFORM CTRANSLATED FROM D.S.D.’S VERSION IN 


NUMSAMPLS = NUMBER OF SEQUENCIES AND SAMPLING INTERVALS 
IN WALSH BASIS CINTEGER) 


DATARRAY = INPUT/OUTPUT ARRAY CREAL) 


68 


o”00 


no 
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10 
20 


30 


10 


29 


30 


ENABLESCALE = TRANSFORM SCALING ENABLE SWITCH 


INTEGER NUMSAMPLS 
REAL DATARRAY (1) 
LOGICAL ENABLESCALE 


CALL BITREVSORT CNUMSAMPLS, DATARRAY ) 
CALL DECINSEG CNUMSAMPLS, DATARRAY > 


IF CENABLESCALE) CALL SCALE (CNUMSAMPLS, DATARRAY ) 
RETURN 
END 


SUBROUTINE BITREVSORT CNUMSAMPLS, DATARRAY ) 
BINARY BIT -— REVERSAL SORT 


INTEGER NUMSAMPLS, NUMSD2, DIRECT, REVERSED, DATAD, DATAR, OF FSET 
REAL DATARRAY(1),SCRATCH 


NUMSD2 = NUMSAMPLS / 2 
REVERSED ® O 
DO 30 DIRECT = 0, NUMSAMPLS - 1 
IF (DIRECT .GE. REVERSED) GO TO 10 
DATAD = 1 + DIRECT 
DATAR = 1 + REVERSED 
SCRATCH = DATARRAY (DATAD) 
DATARRAY CDATAD) = DATARRAY (CDATAR) 
DATARRAY CDATAR) = SCRATCH 
OFFSET = NUMSD2 
IF CCOFFSET .GT. REVERSED) .OR. (COFFSET .LT. 2)) GO TO 30 
REVERSED = REVERSED - OFFSET 
OFFSET = OFFSET / 2 
GO TO 20 
REVERSED = REVERSED + OFFSET 
RETURN 
END 


SUBROUTINE DECINSE@ CNUMSAMPLS, DATARRAY ) 
DECIMATION IN SEQUENCE 


INTEGER NUSAMPLS, LOGZNUMSAMPLS, MSE, RLOCKSIZE, STAGE 
INTEGER NODE1START, NODE2ZSTART,NODE1, NODE2 

LOGICAL TOGGLE 

REAL DATARRAY(1),SCRATCH 


LOGZNUMSAMPLS = O 

MSB = NUMSAMPLS 

MSE = MSE / 2 

IF (MSE .E@. 90) GO TO 20 
LOG2ZNUMSAMPLS = LOG2NUMSAMPLS + 1 
GO TO 10 

BLOCKSIZE = NUMSAMPLS / 2 

DO 80 STAGE = 1, LOG2NUMSAMPLS 
TOGGLE = .FALSE. 
NODEISTART = O 


IF CNODEI1START .GE. NUMSAMPLS) GO TO 80 
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NODE2ZSTART = NODEISTART + BLOCKSIZE 
NODE2 = NODE2START 
IF (.NOT. TOGGLE) GO TO 50 
DO 40 NODE1 = NODEISTART,NODE2START - 1 
SCRATCH = DATARRAYCNODE1 + 1) 
DAT ARRAY (NODE 1+1)=DATARRAY (NODE 1+1)-DATARRAY (NODE 2+1) 
DATARRAY (NODE2+1) = DATARRAY (NODE2+1) + SCRATCH 
NODE2 = NODE2 + 1 
GO TO 70 
DO 60 NODE1 = NODEISTART,NODE2START - 1 
SCRATCH = DATARRAY (NODE1+1) 
DATARRAY (NODE 1+1) =DATARRAY (NODE 1+1)+DATARRAY (NODE2+1 ) 
DATARRAY (NODE2+1) = SCRATCH - DATARRAY (NODE2+1) 
NODEZ = NODE2 + 1 
TOGGLE = .NOT. TOGGLE 
NODE1START = NODE2 
GO TO 30 
BLOCKSIZE = BLOCKSIZE / 2 
RETURN 
END 


SUBROUTINE SCALE (CNUMSAMPLS, DATARRAY ) 
WALSH TRANSFORM SCALING 


INTEGER NUMSAMPLS, I 
REAL DATARRAY (1),SCALEFACTOR 


SCALEFACTOR = 1.0 /NUMSAMPLS 
DO 10 I = 1, NUMSAMPLS 
DATARRAY(I) = DATARRAY(I) # SCALEFACTOR 
CONTINUE 
RETURN 
CLOSE (11) 
END 
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Figure B-7 Diffraction Pattern for Sequency Six 
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Figure B-18 Diffraction Pattern for Sequency 20 
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Figure B-21 Diffraction Pattern for Sequency 23 
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Figure C-1 Cutter512V 


107 





Figure C-2 Cutter256V 


108 


ES 

SECIS 
ce 
yee 


| 


N 





Figure C-3 Cutter128Vv 


109 


a : SSN ES 
FSAI . ) 


< ~ . es ‘ ~ : : : : ‘ 
Ak SASS Sk : WX ES 


SAN 
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Figure E-1 Cutter512 
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Figure F-1 Lhouse512 
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Figure F-3 Lhouse128 
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Figure I-1 Emblem512 
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APPENDIX J 


lee ede i ee SS RAR SSARAAAAAARAERAEEREAERARKAKKKAKKES / 
/* C512: Conversion of 512 X 512 DT-IRIS Images into Simple Byte Format */ 
*x 


/* 910614DSD Microsoft C5.1 
[RRRAARAARRRKRKERERRRKREKKEKEREERERREEREREEEEEEKEEKEEEREEREEREREKEREEEEREEEEEEEEKE / 


#include <stdio.h> 
#include <stdlib.h> 


main (void) 
{ unsigned int a, row, column; 
char instreamname(40], outstreamname[40}; 
x x 
Beton ere outstreanm; 
for (a = 1; a <= 79; a ++) 
mech (¢#?! ) 3 
printf C™\M\E") 3 
for (a = 1; a <= 17; a ++) 
wecn (% *); 
printf ("T512: 512 X 512 DT-IRIS Image Conversion\n\r") ; 
= ° = ; ++ 
Ber ap aie <= 79; a ) 
printf ("\n\n\r") ; 
printf (ea ee DT-IRIS input file ? "); 
scanf ("%s" nstreamname) ; 
rf ((instream = fopen (instreamnane, "rb")) == NULL) 
{ fprintf (stderr, "Can’t open input file!\n") ; 
exit (0); } 
printf CS \n\n\r") ; 
prantt Nena 3 translated output file ? "); 
scanf ("%s outstreamname) ; 
if ((outstream = fopen (outstreamnane, “wb")) == NULL) 
{ exit (0); }. "Can’t open output file! \n"); 
exit (0 } 
peantr ("\n\n\r") ; 
for (a = 1; a <= 512; a ++) 
fgetc (instream); 
for (row = 1; row <= 512; row ++) 
{ for (column = 1; column <= 512; column ++) 
{ a = fgetc (instrean) ; 
irem(aiae, outetreaniis} 
utstreanm) ; 
fclose (instream) ; 
fclose (outstrean) ; 
printf ("Output file contains 512 lines of 512 single-byte unsigned "); 
eee ase 
; 


Ly 


APPENDIX K 


[ RRRREREREREEEEEEEREEEREREREREEAERER ERSTE ee SS Oe 
/* D512: VGA Display for 512 X 512 Single-byte Greyscale Images 


/* 910614DSD Microsoft C5.1 
| Bh te te he te te te te te hh RE RER RRR EEE EERE REEREEEREEEEREREREREREKEEREKEKEKAKEKKEKE 


#include 
finclude 
finclude 
finclude 
#include 


<stdio.h> 
<stdlib.h> 
<conio.h> 
<graph.h> 
<math.h> 


main (void) 
{ int basis, a, amin, amax, row, column, rowoffset, xstart, y, x, color; 
long int tints[16]; 
char instreamname[40); 
FILE * instream; 
struct videoconfig vc; 
basis = 512; 
tints[(0] = 0x00000000; 


tints[1] = 0x00040404; 
tints{2] = 0x00080808; 
tints[{3] = 0x000c0c0c; 
tints[4] = 0x00101010; 
tints[(5]) = 0x00141414; 
tints[6] = 0x00181818; 
tints(7) = Ox00lclcilc; 
tints(8] = 0x00202020; 
tints[(9] = 0x00242424; 
tints[10) = 0x00282828; 
tints[(11]} = Ox002c2c2c; 
tints(12] = 0x00303030; 
tints[(13]) = 0x00343434; 
tints[(14] = 0x00383838; 
tints[(15] = 0x003c3c3c; 
system ("CLS"); 
for (a = 1; a <= 79; a ++) 
putch (’*'); 
prints ¢“\n\5")s 
for (a = 1; a <= 20; a ++) 
putch (’ °); 
printf ("D512: 512 X 512 Greyscale Image Display\n\r") ; 
for (a = 1; a <= 79; a ++) 
utch Me) 


printf ("\n\n\r") ; 
printf ("Name of input file ? "); 
scanf ("%&s", instreamname) ; 
if ((instream = fopen (instreamname, "rb")) == NULL) 
{ fprintf (stderr, "Can’t open input data file! \n"); 


exit (0); 


} 
printf ("\n\n\r"); 


amin 


= amax = 0; 


for (row = 1; row <= basis; row ++) 
{ for (column = 1; column <= basis; column ++) 


{ a = fgetc (instream); 
if (a < amin) 
amin = a; 
else if (a > amax) 
amax = a; 


} 
fgetc (instream); } 


rewind (instream) ; 
setvideomode ( _VRES16COLOR) ; 
_remapallpalette (tints) ; 
_getvideoconfig (& ve 
~Clearscreen ( “GCLEARSCREEN) ; 
rowoffset = (basis - vc.numypixels) / 2; 
for (row = 1; row <= rowoffset; row ++) 


158 


for (column = 1; column <= basis + 1; column ++) 
a = fgetc (instream); 
xstart = (vc.numxpixels - basis) / 2; 
y = 0; 
for (row = 1; row <= vc.numypixels; row ++) 
{ x = xstart; 
for (column = 1; column <= basis; column ++) 
{ a = fgetc (instream) ; 


color = tae ((ve.numcolors - 1) * (a - amin) / (amax - amin) + 


+ 0.5); 
setcolor (color); 
_setpixel (x ¥) +} 

++ xX; 
fgetc (instream) ; 
a Vir 
fclose (instream) ; 
while (! kbhit ()) 


_clearscreen ( GCLEARSCREEN) ; 


_setvideomode (_DEFAULTMODE) ; 
return; } 
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APPENDIX L 


LT B.J. MUSSELMAN 


THESIS PROGRAM —- IMAGE.FOR 


THIS FROISRAM IS CONFIGURED TO PERFORM A 2-DIMENSIONAL WALSH 
TRANSFORM ON AN INPUT IMAGE FILE. THE INVERSE TRANSFORM IS 
THEN TAKEN AND THE TRANSFORMED DATA IS WRITTEN TO SEPARATE 


OUTPUT IMAGE FILE. 


INTEGER*2 NUMROWS, NUMCOLS, MAXSEQ, CUTOFF, SEQUENCY, FILTER(513) 
INTEGER*#2 ROW, COL 

INTEGER#1 BDAT(512,513) 

INTEGER#4 W(513) 


NUMF'OWS Siz 
NUMIZOLS 213 
MAXSEQ = NUMROWS - 1 


OPEN(11, FILE = ’B:\EMBLEM’, FORM=’ BINARY’, STATUS = "OLD’) 
OPEN( 12, FILE=’C:\STUDENTS\ BRIAN\ IMAGE’ , FORM=’ BINARY’, STATUS=’ OLD’ ) 


WRITEC#,#) "ELIMINATE SEQUENCY COEFFICIENTS BEYOND - ° 
READC#,#) CUTOFF 


SEQUENCY FILTER DEFINED IN THIS SECTION 


DO 2 SEQ = 0, MAXSEQ 
FILTER(SEQ + 1) = 1 
IF (SEQ .GT. CUTOFF) FILTER(SEQ + 1) = 0 
CONTINUE 


READ IMAGE FILE, ROW-RBY-ROW INTO ARRAY BDAT 


DO 4 ROW = 1, NUMROWS 
DO 3 COL = 1, NUMCOLS 
READ(11) BDAT CROW, COL) 
CONTINUE 
CONT INUE 


THIS SECTION PERFORMS VERTICAL SEQUENCY FILTERING 


DO 8 COL = 1, NUMCOLS - 1 
WEI TEC#,#) COL 
DO S ROW = 1, NUMROWS 
WCROW) = BDAT CROW, COL) 
CONTINUE 


CALL FWT (NUMROWS, W, .FALSE.) 
DO 6 SEQ = 0, MAXSEOQ 
WC(SEQO + 19 = WCSEQO + 1) #* FILTER CSEQ + 1° 
CONT INUE 
CALL FWT CNUMROWS, W, . TRUE.) 
DO 7 FOW = 1, NUMRFOWS 
BDATCROW,COL) = WCROW) 


CONT INUE 
lONT INUE 
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mY 


15 


16 


17 
18 


19 
1 


THIS SECTION PERFORMS HORIZONTAL SEQUENCY FILTERING 


DO 18 ROW = 1, NUMROWS 
WRITEC#,#) ROW 


DO 15 COL = 1, NUMCOLS 
WCCOL) = BDATCROW, COL) 
CONTINUE 


CALL FWT CCNUMCOLS - 1), W, .FALSE.) 


DO 16 SEQ = O, MAXSEQ 
WCSE@ + 1) = WCSE@ + 1) # FILTERCSEQ + 1) 
CONT INUE 


CALL FWT (C(NUMCOLS - 1), W, . TRUE.» 
DO 17 COL = 1, NUMCOLS - 1 
BDAT(ROW,COL)? = W(COL) 
CONTINUE 
CONTINUE 
WRITE IMAGE FILE, ROW-BY-ROW 
DO 11 ROW = 1, NUMROWS 
DO 19 COL = 1, NUMCOLS 
WRITE(12) BDAT CROW, COLD 
CONTINUE 
CONTINUE 
CLOSE (11) 
CLOSE (12) 
STOP 
END 


SUBROUTINE FWT CNUMSAMPLS, DATARRAY, ENABLESCALE ) 
FAST WALSH TRANSFORM CTRANSLATED FROM D.S.D.°S VERSION IN "C") 


NUMSAMPLS = NUMBER OF SEQUENCIES AND SAMPLING INTERVALS 
IN WALSH BASIS CINTEGER) 


DATARRAY = INPUT/OUTPUT ARRAY CREAL) 
ENABLESCALE = TRANSFORM SCALING ENABLE SWITCH 
INTEGER*#2 NUMSAMPLS 

INTEGER #4 DATARRAY (1) 

LOGICAL ENABLESCALE 

CALL BITREVSORT CNUMSAMPLS, DATARRAY ) 

CALL DECINSECQ CNUMSAMPLS, DATARFAY) 

IF CENABLESCALE) CALL SCALE (NUMSAMFLS, DATAFFAY > 
RETURN 

END 

SUBROUTINE BITREVSORT CNUMSAMPLS, DATARRAY ) 


BINARY BIT — REVERSAL SORT 
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*O * 


10 


30 


10 


20 


30 


40 


30 


60 
7Q 


INTEGER NUMSD2, DIRECT, REVERSED, DATAD, DATAR, OFFSET 
INTEGER#2 NUMSAMPLS 

REAL SCRATCH 

INTEGER#4 DATARRAY (1) 


NUMSD2 = NUMSAMPLS / 2 
REVERSED = O 
DO 30 DIRECT = O, NUMSAMPLS - 1 
IF (DIRECT .GE. REVERSED) GO TO 10 
DATAD = 1 + DIRECT 
DATAR = 1 + REVERSED 
SCRATCH = DATARRAY (CDATAD) 
DATARRAY (DATAD?) = DATARRAY (DATAR? 
DATARRAYCDATAR) = SCRATCH 
OFFSET = NUMSD2 
IF CCOFFSET .GT. REVERSED) .OF. COFFSET .LT. 29) GO TO 30 
REVERSED = REVERSED - OFFSET 
OFFSET = OFFSET / 2 
GO TO 20 
REVERSED = REVERSED + OFFSET 
RETURN 
END 


SUBROUTINE DECINSEQ (CNUMSAMPLS, DATARRAY ) 
DECIMATION IN SEQUENCE 


INTEGER LOGZ2NUMSAMPLS, MSB, BLOCKSIZE, STAGE 
INTEGER NODE1ISTART, NODE2ZSTART, NODE1,NODE2 
INTEGER#2 NUMSAMPLS 

INTEGER#4 DATARRAY (1) 

LOGICAL TOGGLE 

REAL SCRATCH 


LOGZNUMSAMPLS = 0 
MSB = NUMSAMPLS 
MSB = MSE / 2 
IF (MSB .EQ@. O) GO TO 20 
LOG2NUMSAMPLS = LOGZNUMSAMPLS + 1 
GO TO 10 
BLOCKSIZE = NUMSAMPLS / 2 
DO 80 STAGE = 1, LOG2NUMSAMPLS 
TOGGLE = .FALSE. 
NODEISTART = O 
IF (NODEISTART .GE. NUMSAMPLS) GO TO 80 
NODE2START = NODEISTART + BLOCKSIZE 
NODE2 = NODE2START 
IF ¢€.NOT. TOGGLE) GO TO So 
DO 40 NODE1 = NODEISTART,NODEZSTART - 1 
SCRATCH = DATARRAYCNODE1 + 1) 
DATARRAY (NODE 1+1)2=DATARRAY (NODE 1+1)-—-DATARRAY (NODE2+1) 
DATARRAY (NODEZ+1) = DATARRAY (CNODE2+1)2 + SCRATCH 
NODEZ = NODE2 + 1 
GO TO 70 
DO 60 NODE1 = NODE1START,NODE2START - 1 
SCRATCH = DATARRAY (NODE1+1) 
DATARRAY (NODE 1+1)=DATARRAY (NODE 1+1)+DATARRAY (NODE2+1) 
DATARRAY (NODE2+1) = SCRATCH —- DATARRAY (NODE2+1) 
NODE2 = NODE2 + 1 
TOGGLE = .NOT. TOGGLE 
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80 


10 


NODEISTART = NODE2 

GO TO 30 
BLOCKSIZE = BLOCKSIZE / 2 
RETURN 
END 


SUBROUTINE SCALE (CNUMSAMPLS, DATARRAY ) 
WALSH TRANSFORM SCALING 


INTEGER 1 

INTEGER#2 NUMSAMPLS 
INTEGER#4 DATARRAY (1) 
REAL SCALEFACTOR 


SCALEFACTOR = 1.0 /NUMSAMPLS 
DO 10 I = 1, NUMSAMPLS 
DATARRAY CI) = DATARRAYCI) *# SCALEFACTOR 
CONT INUE 
RETURN 
END 
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